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Radiative corrections to b — > crv T 
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Abstract 

Analytical calculation is presented of the QCD radiative corrections to the rate of the 
process b — > ctv t and to the r lepton longitudinal polarization in tv t rest frame. The 
results are given in the form of one dimensional infrared finite integrals over the invari- 
ant mass of the leptons. We argue that this form may be optimal for phenomelogical 
applications due to a possible breakdown of semilocal hadron-parton duality in decays 
of heavy flavours. 



The semileptonic decay rate of B mesons is one of the key ingredients in the deter- 
mination of weak mixing angles. Transitions of the b quark to the charmed as well as 
to the up quark have been analysed in great detail, exploiting either the inclusive decay 
rate or exclusive channels. The analysis of the inclusive rate is, however, affected by the 
uncertainty in the b mass and by bound state corrections. This problem is only partly 
circumvented by fixing m& — m c through the difference of bottom and charmed meson 
masses and by relating the bound state effects to phenomenological constants that can 
be determined from other observables in the context of the heavy quark effective theory 
(HQET). The decay rate is, furthermore, affected by perturbative QCD corrections, 
which have been calculated analytically up to order a s for arbitrary b and c masses 
and massless leptons [||, |], [|. Numerical results for perturbative QCD corrections to 
the partial decay rate b — > rvX have been recently obtained in Ref.[|J. The bound 
state corrections to this decay chanel up to order 1/m 2 are also known ||, ||, U in the 
context of the HQET. The comparison between decays into light (/x, e) and heavy (r) 
leptons may furthermore help to test the theoretical approach and allow to fix some of 
the free parameters. Including b — > u transitions, four kinematically different leptonic 
decay modes are thus available for the comparison. 

In this paper analytical results for the decay rate b — > crv T are presented in the 
form of an one dimensional integral over the invariant mass squared w 2 of the leptonic 
system. This formulation allows, at least in principle, the separation of the region 
of relatively small w 2 , where the inclusive parton model description based on local 
parton-hadron duality should work, from the region of large w 2 where only one or few 
resonances are produced and the duality between the parton and hadron description 
may be doubtful^. In the region of large w 2 i.e. close to the Shifman-Voloshin limit || 
w 2 = w^ aa . the theoretical description of lepton spectra based on summation over 
exclusive channels is particularly simple and reliable. The HQET approach for exclusive 



decays 0, [10| on the other hand becomes quite involved if not impractical in the region 
of small w 2 which is dominated by multiparticle final states. 

The calculation of the lowest order rate as well as corrections can be related in a 
straightforward way to the corresponding calculations for the decay into a virtual W 
boson with the subsequent integration over the mass of the Iv system. The differential 
decay rate is proportional to 

H af3 £ Q/3 dPS (b ^ctv) 

Tt a p depends on quark and gluon fields and 

£ a f3(r;is) ~ ^[w r 7 a (l - J 5 )v u ] [u T 7/3(l-75K ] f 

s 
~ V a Tp + T a Up - V ■ rg a p - \£ a p lS V t T £ (1) 



^n a recent preprint M the breakdown of the local parton-hadron duality has been invoked as the 
origin of problems with the semileptonic decay rate of D mesons in the framework of HQET. Let us 
remark that for large w 2 the kinetic energy of the hadronic system in B decays can be similar to that 
in D decays. Thus the semileptonic branching ratios for b decays may be also affected for the effective 
mass of the hadronic system close to the resonance region. 



where r a and v a are the four- momenta of r and v T . The phase space for the decay 
of b into ctv is, in the standard way, decomposed into a sequence of two-particle final 
states 

dPS(6 -► crv) ~ dw 2 dPS(6 -»■ cw) dPS(w -► rz/) (2) 

where w a = r a + v a . Then, it is straightforward to show that 

|dPS(w - ru)C aP ~ A (m 2 /w 2 ) T$ + B (m 2 /w 2 ) 7$ = 

( 1 - m 2 /w 2 f [ ( 1 + 2m 2 /w 2 ) w qW/3 - ( w 2 + \m\ ) ^] (3) 

where 

rp(0) 

T <$ = w « w /3 - w 2 ^ (4) 

It follows that the decay rate can be split accordingly into two incoherent terms which 
are related to weak decays of a heavy quark Q into another quark q and a real spin one 
or a spin zero boson. The relative weight of spin one versus spin zero contributions is 
governed by their respective spectral functions and it can be derived from eq.(^|). For 
massless leptons .4.(0) = and 13(0) = 1, and therefore only the spin one (transversal) 
component (~ T^J) contributes. The result is given in 0. For fixed w 2 this contribu- 
tion to the rate can be obtained from the formula for t — > bW, the top quark decay 
into b quark and a real W boson f||fj Multiplying this formula by B (m 2 /w 2 ) one 
obtains the contribution of the spin one component for m T ^ 0. The other (longitudi- 
nal) contribution (~ T^J) can be in an analogous way related to the (yet unobserved) 
decay t — > bH + where H + denotes a charged Higgs boson. Let Q and q denote the 
four-momenta of the quarks in the initial and in the final state. For a two-body decay 
mode the momentum of the W boson is W = Q — q. In Born approximation 

W,u(q)^(l- l5 )u(Q) = u(q)[(Q-q)(l- l5 )]u(Q) 

= u(q)[(M-m) + (M + m)j 5 ]u(Q) (5) 

where the equations of motion 

(Q - M) u(Q) = and u(q) (q - m) = (6) 

have been used. The last line in (||) can be interpreted as the amplitude of the decay 
Q —> qH where if is a spin zero boson whose coupling to the weak quark current is 
given by 

g = (M -m) + (M - 772)75 (7) 



2 Note that the rates of up and down type quark decays into their respective isospin partners are 
of course identical, in contrast to the shapes of the lepton spectra. 



Although not applicable for individual Feynman diagrams, the same relation holds 
true for the longitudinal contribution (~ W^) to the decay amplitude when 0(a s ) 



QCD corrections are included [II]. Therefore this contribution to the rate b — > crv T 



can be extracted from a formula describing t — > bH + which has been given in [12|; cf. 

Model I therein. 

Let us define now the following dimensionless quantities 





p = m 2 Jm 2 b 


1] = m^/m 2 . 


and the kinematic functions 
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where w" = r a + v a and w 2 denotes the effective mass of rv T . The partial rate of the 
decay b — > ctv t is given by 



CTV T 



(i-Vp) 2 dr 



i\t 



dt 



(9) 



with the differential rate 
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(10) 
(11) 



Fo(t) 



4p 3 
4p 3 



(l-p) 2 + £(l + p)-2£ 2 
(l-p) 2 -t(l + p) 



A{$ + A 2 Y W + A 3 Y P + Am lnp + A bP3 
B^ + B 2 Y W + B 3 Y P + B 4 p 3 lnp + B 5 p 3 



(12) 

(13) 
(14) 
(15) 



^ = 81n(2p 3 ) - 2 hit + [2Li 2 (w_) - 2Li 2 (w + ) + 4Li 2 (2p 3 /p + ) 

— 4Ypln(2p 3 /p + ) — lnp_ \nw + + lnp + lnw_] 2p /p 3 



(16) 
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p) l+£-4£ 2 -p(2-£) + p 2 
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B, 
B 5 



3 + 6t - 21£ 2 + 12t 3 - p(l + 12* + 5£ 2 ) + p 2 (ll + 2t) - p 3 
1 + 3t - 4£ 2 - p(4 - t) + 3p 2 " 
5 + 9£ - 6£ 2 - p(22 - 9t) + 5p 2 " 

-8(i-p)[(i-p) 2 -t(i+py 

-4(1 - p) 4 /t - 2(-l + 3p + 15p 2 - 5p 3 ) + 8(1 + p)t - 6(1 + p)t 2 
-4(1 - pf/t - 2(1 - p)(l - lip) + 6(1 + 3p)t 
-6(l-3p)(3-p) + 18t(l + p) 



(17) 



The integral in eq.(|9|) can be easily performed for the Born contribution. It reads: 



r (b 



CTV T 



v hc { 24 [ v 2 (i - p 2 ) y w + p 2 (i - v 2 ) y p 



+ 



2V 3 (l - 7rj - 7r] 2 + rf - 7p + 12r/p - 7r/ 2 p - 7p 2 - 7r/p 2 + p 3 ) } 



(19) 



where 



^3 =Pz{v), 



y P = ym 



and 



y w = y w ( v ). 



In principle the first order QCD correction can be also expressed in terms of polylog- 
arithms. In particular for 77 = the formula (12) of [0] is obtained. However, the 
complete result is lenghty. From the practical point of view it is much simpler to eval- 
uate the integral in eq.([|) numerically. Moreover, as it has been explained, for t close 
to the Shifman-Voloshin limit the exclusive description is preferable, thus, for applica- 
tions, it may be better to perform this integral only over a part of the available phase 
space. In recent articles |15|, [RJ the size of 0(a 2 ) corrections has been estimated using 
the scheme of Brodsky, Lepage and Mackenzie |17| for fixing the scale p of a s . It has 
turned out that this scale is rather low. This suggests that next-to-leading QCD cor- 
rections are large. The most serious problems arise in the region close to the no-recoil 
t = t max point. It is well known, cf. |L8], |H|, that at the boundaries of the available 
phase space logarithmic divergences may appear as remnants of infrared divergences. 
This phenomenon arises when real radiation becomes supressed relative to virtual cor- 
rections. This is exactly what happens for t in the neighbour of t max and once again 
one is led to the conclusion that the inclusive parton-like description may break down 
there. 

In the massless limit p — >• which corresponds to b —>■ utv t transition the functions in 
eq.(|15D simplify considerably 



Fo{t) 


= 2(1 - t) 2 (l + 2t) 


(20) 


m) 


= 2(1 -t) 2 


(21) 


r x {t) 


= F (t) [§7r 2 + 4Li 2 (t) + 21ntln(l-t)] - (1 - £)(5 + 9i - 6t 2 ) 






+ 4t(l - t - 2£ 2 ) In t + 2(1 - t) 2 (5 + At) ln(l - t) 


(22) 



F»{t) = f*(t) [|tt 2 + 4Li 2 (t) - | + ln(l-t)(2 hi£-f + 5)] +4(1 - t)t hit 

(23) 

After integration over t one derives the following expression for the total partial rate 
of b — > urv T 

r(6 -> urz/ r ) = 1 - 877 + 8?7 3 - if - 12r] 2 \n7] 

— — ((-1 + n)(75 - 539?7 - 476r] 2 + 18r] 3 )/12 

+ (3 - 24r/ - 36r/ 2 + 16r/ 3 - 2?7 4 )7r 2 /3 

+72t/ 2 [C(3) - U 3 (t])} + 2(1 -877 + 36r/ 2 + 16r? 3 - 2r ] i )U 2 (r ] ) 

+ (1 - r/ 2 )(31 - 320?7 + 31?7 2 ) ln(l - 77) /6 



+ 



2r] + 15r/ 2 - 94r? 3 /3 + 31r/ 4 /6 - 8r/V + 24r7 2 Li 2 (r7) 

+2(1 - rf)(l -8rj + r/ 2 ) ln(l - 77)] ln(r/)} (24) 



G 2 m5|K.J 2 



6m 



(25) 
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which agrees with eqs.(2.21) and (3.6) in a recent preprint 

It has been argued |l | that the r polarization in b — > ctv t is particularly sensitive 
to deviations from the Standard Model. The longitudinal component of r polarization 
can either be measured with reference to its direction of flight in the b rest frame or, 
alternatively, relative to its direction of flight in the rv rest frame. To evaluate an- 
alytically QCD corrections to the longitudinal r polarization in the b rest frame is a 
demanding task. Previous experience from a similar calculation for the decay of polar- 
ized top quarks indicates that these corrections are typically quite small through most 
of the kinematical range. To substantiate this claim, the longitudinal r polarization in 
the rv rest frame is evaluated including QCD corrections. 
For V-A weak current the leptonic tensor 

C af 3(r,s;u) ~ [u T j a (l -75)^] [u T ^p{l - ^)v v } ] 

~ v a T p + T a pp - v ■ Tg aP - ie al3y s^'T 8 (26) 

where 

T a = \(T a -m T s a ) (27) 

and s a is the r polarization fourvector. The helicity states in the rv rest frame are 
obtained for 

m T s a = ±(r a - ^v a ) (28) 

It is evident that the net polarization can be calculated in the same way as the total rate, 
decomposing again the spin dependent term into longitudinal and transversal parts. 



For the positive helicity of r one derives the following expression for the differential 
rate 

£ ^r = I ( l ~ f) {*m + 3 ^w - w [*«> + 3 ^ (i) '} (29) 

For the negative helicity one has 

dr(-) dr drw 

(30) 



dt dt dt 

and the net helicity in the tv t rest frame is equal to 

P--1 + Zi2 df (oi\ 

h At ut 

In the subsequent discussion the mass difference m& — m c will be fixed through the 
HQET relation 

m b -m c = [(3m B * + m B ) - (3m D » + m jD )]/4 + . . . 

The most important corrections to the above relation arise from the kinetic energy of 
heavy quarks in B and D mesons. On physical grounds one expects this contribution 
to increase the difference between mj and m c . We take mt, — m c = 3.4 GeV in numerical 
calculations and mj is varied between 4.5 and 5.0 GeV. In order to estimate effects of 
QCD corrections on measurable quantities we neglects all the problems with the parton- 
hadron duality and the scale ambiguity. The following results have been obtained: 

r(6 -»• crv T ) = T (b -► ctv t ) [ 1 + (-0.450 ± 0.016) a s ] 
Y(b^cev e ) = T (b^ceu e ) [1 + (-0.545 ± 0.025) a s ] 
R = BR(b -> rX)/BR(b -> eX) = R [ 1 + (0.094 t 0.009) a s ] 
P = -(0.293 =p 0.002) [ 1 + (0.140 T 0.015) a s ] 

< t > (b -► tX) = (0.34 ± 0.03) [ 1 + (0.016 ± 0.003) a s ] 

< t > (b -»• eX) = (0.20 ± 0.02) [ 1 + (0.035 ± 0.007) a s ] 

where < t >(...) denote the average values oft for the corresponding decay chanels. It 
is evident that the 0(a s ) corrections practically cancel in the ratio R of the branching 
ratios as well as in the result for the polarization P. The moments < t > (b — > rX) 
and < t > (b — > eX) are also insensitive to a s corrections. On the other hand all these 
quantities are sensitive to the quark masses and therefore they may be used for fixing 
nib and m c . 
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